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Abstract 

In this paper we study the rate of convergence of a symmetrized version of the Milstein scheme applied to the 
solution of the one dimensional SDE 

X t = x 0 + f b(X a )ds + [ a\X 3 \ a dW s , x 0 > 0, <7 > 0, ct € [§, 1). 

Jo Jo 

Assuming 6(0) /a 2 big enough, and 6 smooth, we prove a strong rate of convergence of order one, recovering the 
classical result of Milstein for SDEs with smooth diffusion coefficient. In contrast with other recent results, our proof 
does not relies on Lamperti transformation, and it can be applied to a wide class of drift functions. On the downside, 
our hypothesis on the critical parameter value 6(0)/cr 2 is more restrictive than others available in the literature. Some 
numerical experiments and comparison with various other schemes complement our theoretical analysis that also 
applies for the simple projected Milstein scheme with same convergence rate. 


1 Introduction and main result 

The Milstein scheme was introduced by Milstein in [16] for one dimensional Stochastic Differential Equations (SDEs) 
having smooth diffusion coefficient. Introducing an appropriated correction term, this scheme has better convergence 
rate for the strong error than the classical Euler-Maruyama scheme. Typically, when the drift and diffusion coefficient 
of one dimensional SDE are twice continuously differentiable with bounded derivatives, the Milstein scheme is of 
order one for strong error (see eg. Talay [20]) instead of one-half for the Euler-Maruyama scheme. This well-know 
fact produces remarks on blogs, internet forums, and software packages that sometimes recommend to use the Milstein 
scheme for constant elasticity of variance (CEV) models in finance, or its extension with stochastic volatility as SABR 
model, (see e.g Delbaen and Shirakawa [ 8 ] and Lions and Musiela [15] for a discussion on the (weak) existence of 
such models); CEV are popular stochastic volatility models of the form 


dX t = fj,X t dt + aX?dW t 

with 0 < 7 < 1. But the interesting fact in this story is that the rate of convergence of the Milstein scheme, for such 
family of processes with 0 < 7 < 1 is not yet well studied, to the best of our knowledge. 
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In this paper we establish a rate of convergence result for a symmetrized version of the Milstein scheme applied to 
the solution of the one dimensional SDE 

Xf = Xq + f b(X s )ds + 

Jo 


a\X s \ a dW s 


( 1 . 1 ) 


where Xq > 0, a > 0 and | < a < 1. Of course Equation (1.1) does not satisfies the hypothesis to apply the 
classical result of Milstein [16], In particular, the diffusion coefficient is only Holder continuous whereas the classical 
hypothesis is to have a C 2 diffusion coefficient. 

The main picture of our convergence rate result is that Milstein scheme stays of order one in the case of Equation 
(1.1), but some attention must be paid to the values of 6(0), a and a. 

There exist in the literature other strategies for the discretization of the solution to (1.1). There are some results 
based on the Lamperti transformation of the equation, for example, by Alfonsi [1, 2], and by Chassagneux, Jacquier 
and Mihaylov [6], And also, there some are results where the equation (1.1) is discretized directly, as in Berkaoui, 
Bossy and Diop [3] or in Kahl and Jackel [12]. In the numerical experiments section, we compare the symmetrized 
Milstein scheme with a selection of schemes proposed in the aforementioned references. We also experiment the 
symmetrized Milstein scheme in a multilevel Monte Carlo application and we compare with other schemes. 


In the whole paper, we work under the following basis-hypothesis: 

Hypothesis 1.1. The power parameter a in the diffusion coefficient of Equation (1.1) belongs to [2, 1). The drift 
coefficient b is Lipschitz with constant K > 0, and is such that 6(0) > 0. 

Hypothesis 1.1 is a classical assumption to ensure a unique strong solution valued in R + . We assume it in all the 
forthcoming results of the paper, without recall it explicitly. To state the convergence result (see Theorem 1.6), another 
Hypothesis 1.5 will be added and discussed, that in particular constrains the values a, 6(0) and a. 


1.1 The symmetrized Milstein scheme 

To complete our task we follow the ideas of Berkaoui, Bossy and Diop in [3] who analyze the rate of convergence 
of the strong error for the symmetrized Euler scheme applied to Equation (1.1). Although, whereas they utilize an 
argument of change of time, we consider first a weighted L p (f2)-error for which we prove a convergence result, and 
then we utilize this result to prove the convergence of the actual //'(O)-error. 

We consider Xq > 0, T > 0, and N £ N. We define the constant step size At = T/N and tk = kAt. Over this 
discretization of the interval [0, T] we define the Symmetrized Milstein Scheme (SMS) (X tk , k = 0,..., N) by 


Xt k = 


x 0 , for k = 0, 

2 

X tk _ ± + b(X tk _JAt + crXt k i {Wt k - W tk _ 1 ) + ^-X 2 C! [(W tk - W t fe _J 2 - At] 


for k = 1,..., N. 

In the following, we use the time continuous version of the SMS, ( X t , 0 < t < T) satisfying 


X t = 


x v{t) + b(X v{t) )(t - fff)) + ax; {t) (w t - W v{t) ) + °^~Xl a {t) 1 [(Wt - W v(t) ) 2 - (t - fft))] 


, ( 1 - 2 ) 


where p(t) = sup^.g^ Ar}{^fc : tk < t}. We also introduce the increment process ( Z t , 0 < t < T) defined by 

2 

Zt = X„ (t) + b(X v(t) )(t - ?7(f)) + aX; (t) (W t - W vW ) + [(W t - W v{t) f - (t - 77 (f))] , (1.3) 

so that X t = \Zt\- Thanks to Tanaka’s Formula, the semi-martingale decomposition of X t is given by 


X t = x 0 + / sgn(Z s )b(X v{s) )ds + L° t (X) + I sgn(Z s ) trX“ (s) + acr 2 X^ s) 1 (W S - W v{s) ) 


dW s (1.4) 


where sgn(z) = 1 - 21[ x < 0 ]. 
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Moment upper bound estimations for X and X 

We summarize some facts about the process (X t , 0 < t < T), the proofs of which can be found in Bossy and Diop [5]. 


Lemma 1.2. For any q > 1, there exists a positive constant C depending on q, but also on the parameters 6(0), K, a, 
a and T such that, for any Xq > 0, 


E 


sup Xt q 


L 0 <t<T 


<C{ l + x 2q ). 


When \ < a < 1, for any q > 0, 

When a = \, for any q such that 1 < q < — 1 


sup E [Xf q ] <C{l + x^ q ). 

0 <t<T 


(1.5) 


( 1 . 6 ) 


sup E [X t q ] < Cx 0 q . (1.7) 

0 <t<T 


Lemma 1.3. Let (X t ,0 < t < T) be the solution of (1.1) with J < a < 1. For all p > 0, there exists a positive 
constant C{T , p), increasing in p and T, depending also on b , cr, a and Xq such that 


Eexp 



<C{T,p). 


( 1 - 8 ) 


When a = the inequality (1.8) holds ifb{ 0) > a , 2 and p < ^-( 2b J-^ — l) 2 . 

Notice that the condition 6(0) > cr 2 /2 is also imposed by the Feller test in the case a = ^ for the strict positivity 
of X, that allows to rewrite Equation (1.1) as 


X t = x 0 + f b(X s )ds + [ a\fXsdW s . 
J 0 Jo 


Using the semimartingale representation (1.4), we prove the following Lemma regarding the existence of moments 
of any order for X t . 


Lemma 1.4. For any q > 1, there exists a positive constant C depending on q, but also on the parameters 6(0), K , a, 
a and T such that for any Xq > 0, 


E 


^2 q 

sup A t 
te[o,T] 


<C{l+xl q ). 


The proof of this lemma is based on the Lipschitz property of 6 and classical combination of Ito formula and Young 
Inequality. For the sake of completeness, we give a short proof in the Appendix. 


1.2 Strong rate of convergence 

The main result of this works is the strong convergence at rate one of the SMS X to the exact process X. The 
convergence holds in L p for p > 1. To state it, we add to Hypothesis 1.1 the following. 

For any x in R + , we denote [:/;] the rounded up integer. 

Hypothesis 1.5. 

(i) Let p > 1. To control the L p (Li)-norm of the error, if a > \ we assume 6(0) > 2a(l — a) 2 er 2 . Whereas for a = J 
we assume 6(0) > 3(2[p V 2] + l)er 2 /2. 

(ii) The drift coefficient b is of class C 2 (R), and b" has polynomial growth. 
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We now state our main theorem. To lighten the notation, we consider for a £ (1.1) 

( b a (a) := 6(0) — 2(1 — a) 2 aa 2 , 

\ K(a ) := K + ^(2a- 1)[2(1 - a )]~^^ 


and we extend this definitions to a = | taking limits. So, b a ( 1) = lim i b a {a) = 6(0) — er 2 /4, and K(^) = 

a-> 2 

lim i K(a) = K. Notice that limQ,^! K(a) = K + er 2 /2, and since K(a) is continuous on (|, 1), we have that 

a ->2 

K(a) is bounded. This is especially important in the definition of A max (a) bellow, because tells us that a H>■ A max (ct) 
is strictly positive and bounded on [|, 1). 

Theorem 1.6. Assume Hypotheses 1.1 and 1.5. Define a maximum step size A max (a) as 


A max (cr) 


x 0 

(1 - y/a)b a (a) 


A 


1 


4,aK(a) ’ 
4iT A " 0 ’ 


foraG (1,1) 
for a =\. 


( 1 . 10 ) 


Let ( X t , 0 < t < T) be the process defined on (1.1) and (X t , 0 < t < T) the symmetrized Milstein scheme given in 
(1.2). Then for any p > 1 that allows Hypotheses 1.5, there exists a constant C depending on p, T, 6(0), a, a, K, and 
xo, but not on At, such that for all At < A max (a), 

sup (E[|A t —A t |P] V <CAt. (1.11) 

0 <t<T v J 


About Hypothesis 1.5. Notice that for a , Assumption ( i ) does not depend on p and becomes easier to fulfill 
as a increases. On the other hand, for a = b. Assumption (*) depends on p in a unpleasant manner. However, as we 
will see later in Section 4 (see Table 1), this kind of dependence in p is expected, and similar conditions are asked in 
the literature for other approximation schemes in order to obtain similar rate of convergence results. 

Also, notice that (*) is a sufficient condition: in the numerical experiments we still observe a rate of convergence of 
order one for parameters that do not satisfy it, but we also observe that for parameters such that 6(0) <C cr 2 , although 
the convergence occurs, it does in a sublinear fashion. 

On the other hand, Assumption ( ii) is the classical requirement for the strong convergence of the Milstein scheme. 
As we will see later in the proof of the main theorem, with the help of the Ito formula, this hypothesis let us conclude 
that 

E [\X„ - Z s | 2 *'- 1 (6(X j?(s) ) - 6(A' S ))] < C ( supE [\X U - X u \ 2 ?] + At 2 ?) 

\u<s J 

instead of 

E [|A' S - A s | 2 p- 1 (6(A ?)(s) ) - 6(A' S ))] <c( supE [\X U - X u \ 2p ] + At p ) , 

\u<s J 

which is the classical bound obtained for the Euler-Maruyama scheme under a Lipschitz condition for a drift 6. 

The rest of the paper is organized as follow. In Section 2 we state some preliminary results on the scheme which 
will be building blocks in the proof of Theorem 1.6. Section 3 is devoted to the proof of the convergence rate. The 
main idea is first to introduce a weight process in the -error. Wet get the rate of convergence for this weighted 

error process, and we use this intermediate bound to control the L 2p (fl)-error, from where we finally control the 
//'(O)-error. Also, as a byproduct we obtain the order one convergence of the Projected Milstein Scheme (see 3.3). 
In section 4, we display some numerical experiments to show the effectiveness of the theoretical rate of convergence 
of the scheme, but also to test Hypotheses 1.5-(i) on a set of parameters. In this section we also shows how the 
inclusion of the SMS in a Multilevel Monte Carlo framework could help to optimize the computational time of weak 
approximation of assets valuation. In Section 5 we present the proof of the preliminary results on the scheme. Finally 
we have included in a small appendix a couple of proofs to make this paper self contained. 
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2 Some preliminary results for X 


This short section is devoted to state some results about the behavior of X, their proofs are postponed to Section 5. 
All these results hold under Hypothesis 1.5-(«) which is in fact stronger than what we need here. So, we present the 
next lemmas with their minimal hypotheses (still assuming Hypothesis 1.1). 

Lemma 2.1 (Local error). For any Xq > 0, for any p > 1, there exists a positive constant C, depending on p, T, the 
parameters of the model 6(0), K, a, a, hut not on At such that 

sup E[\X t -X v(t) \ 2p ] <CAt p . 

0 <t<T 


By construction the scheme X is nonnegative, but a key point of the convergence proof resides in the analysis 
of the behavior of X or Z visiting the point 0. The next Lemma shows that although Z t is not always positive, the 
probability of Z t being negative is actually very small under suitable hypotheses. 

Lemma 2.2. For a G [ 2 ,1), ifb( 0) > 2ct(l — a) 2 a 2 , and At < 2 K(a) > ^ len there exists a positive constant 7 , 
depending on the parameters of the model, but not on At, such that 


sup J 

k=0,...,N-l 


inf Z s < 0 ) < exp (— . 

t k <s<t k+ 1 / V At) 


In particular, 


sup P ( Z t < 0) < C exp . 

0 <t<T v 


( 2 . 1 ) 


( 2 . 2 ) 


To prove Lemma 2.2, it is necessary to establish before the following one, which although technical, gives some 
intuition about the difference between the SMS and the Symmetrized Euler scheme presented in [3], 


Lemma 2.3. For a G [|, 1), ifb( 0) > 2a(l — a) 2 a 2 , we set x(a) = > 0. Then for all t G [0, T], and for 

P G (0, 1]> _ 

P [Z t < (1 - p)b a {a)At 1 X v ^ < px[a) ] = 0. 

Roughly speaking, from this lemma we see that when Z v (^) > 0, Z t becomes negative only when 


■ all 


I z t - Z v(t) | > px{a). 

But observe that only the left-hand side of this inequality depends on Af, and its expectation decreases to zero propor¬ 
tionally to v Af, according to Lemma 2.1. 

Now imposing At small enough, we prove an explicit bound for the local time moment of X. 

Lemma 2.4. For a G [|, 1), if 6(0) > 2a(l — a) 2 a 2 and At — 2 K(a) ^ (i (a) ’ ^ len there exist positive 
constants C and 7 > 0 depending on a, 6(0), K, and a but not in At such that 


E(L° t (X) 2 ) < exp 

We end this section with another key preliminary result, which is the convergence rate of order 1 for the corrected 
local error. Although the classical local error is of order 1/2, as stated in Lemma 2.1, the local error seen by the 
diffusion coefficient function, corrected with the Milstein term stays of order 1. 

Lemma 2.5 (Corrected local error process). Letusfixp> 1, and a G [ 7 , 1 ). Fora > assume 6(0) > 2a(l— a) 2 a 2 , 
whereas for a = \, assume 6(0) > 3(2 p + l)cr 2 /2. Then, there exists C > 0, depending on the parameters of the 
model but not in At, such that for all At < A max (a), the Corrected Local Error satisfies 



sup E 
0 <t<T 


aX? - aX a v(t) - aa 2 X% t) \w t - W vW ) 


< CAt 2p . 
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3 Proof of the main Theorem 1.6 


The proof of Theorem 1.6 is built in several steps. First, we work with the L 2p (fl)-norm of the error, for p > 1, then 
at the last step of the proof we go back to the L p (fl)-norm for p > 1. 

In what follows we denote 

S t := X t ~ X t 


and 

so that 


£ t := sgn (Z t ) crX* (t) + aa 2 X^ t) (W t - W v{t) ) 


-aX? 


dS t = (sgn (Z t )b(X v(t) ) - b(X t )) dt + dL° t (X) + Z t dW t . 


Also, to make the notation lighter, we will denote the Corrected Local Error by 


D t {X) := - aX^ t) - aa 2 X^ t) \w t - W v(t) ). 


3.1 The Weighted Error 


Before to prove the main theorem, we establish in the Proposition 3.1 the convergence of a weighted error. For p > 1, 
let us consider (/3 t , 0 < t < T), defined by 


Pt = 2p||6'||oo + 2p(4p - 1) + 


8a 2 p(4p — 1 )a 2 

Xf (1 ~ a) 


and the Weight Process (r*, 0 < t < T) defined by 



(3.1) 


(3.2) 


The Weight Process is adapted, almost surely positive, and bounded by 1. Its paths are non increasing and hence 
has bounded variation, and also satisfies 

dT t = -/3 t r t dt. 

The process (r t ) f>(| can be seen as an integrating factor in the sense of linear first order ODE (see for example 
[19]). When we apply the Ito’s Lemma to T 2 £f p , instead of £ f Vp alone, we can remove a very annoying term that 
appears in the righthand side bound (see the proof of Lemma 3.3). The exponential weight in a //'(f>)-norm is a 
useful tool to obtain a priori bound (as an example, for the existence and uniqueness of the solution of a Backward 
SDE, it is introduced a norm with exponential weight, such that the operator associated to the BSDE is contractive 
under this new norm (see proof of Theorem 1.2 in [17]). In the same way, here we introduce this exponential weight 
to get the following a priori error-bound, which will allow us to prove Theorem 1.6. 

Proposition 3.1 (Weighted Error). Under the hypotheses of Theorem 1.6, for p > 1 and a € |3,1), there exists a 
constant C not depending on At such that for all At < A max (a) 


sup E 

0 <t<T 


r isv 


< CAt ip . 


(3.3) 


Remark 3.2. 

(i) Since we are going to work first with the L 2p (fiY)-norm of the error, when a = \, Hypothesis 1.5-(i) becomes 

m > 3(2[2p V 2] + l)cr 2 = 3(4p + l)g 2 ^ 


in particular 



> 12 p + 2. 


(3.4) 
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(ii) From Lemma 1.2, the process (3 has polynomial moments of any order for a > 4, and when a = g, there exists 
C such that E[f3f] < C,for all 1 < q < 2b(0)/cr 2 — 1. From the previous point in this Remark, it follows that the 
process f3 has moments at least up to order 12 p + 2. 

(Hi) From the definition f3 in (3.1), and due to Lemma 1.3, there exists a constant C such that E[r^ 9 ] < C for all 
q > 0 when a > whereas for a= \, the q-th negative moment of the weight process is finite, as soon as 


2p(4p - l)er 2 q<Y 



Notice that, thanks to point (i) in this Remark, a sufficient condition such that this last inequality holds is 


I6p(4p — 1) q < (12p + 2) 2 . 


We cut the proof of Proposition 3.1 in two technical lemmas. 

Lemma 3.3. Under the hypotheses of Theorem 1.6, for p > 1 and a £ [^, 1), there exists a constant C not depending 
on At such that for all At < A max (a) 

E{r 2 t £? p } < 4p f E [r 2 ^- 1 (b(X v(s) ) - b(X s ))] ds + 4p||& , || 00 f sup E[Tl£ 4 u p ]ds + CAt ip . (3.5) 

JO Jo 0<u<s 

Lemma 3.4. Under the hypotheses of Theorem 1.6, for p > 1 and a £ [^, 1), there exists a constant C not depending 
on At such that for all At < A max (a), and for any s G [0, T] 

|E[r 2 £ , s 4p - 1 (6(A J)(s) ) - 6(X S ))]| < CsupE[r 2 ^] + CAt 4p . (3.6) 

lf<S 


As we will see soon, we prove (3.6) with the help of the Ito’s formula applied to b, and here is where we need b of 
class C 2 required in Hypotheses 1.5 -(m). 


Proof of Proposition 3.1. Thank to Lemmas 3.3 and 3.4, we have 

E[T 2 t £? p ) <C f supE[r 2 £^]d S + CAt 4p . 

Jo u<s 

and since the right-hand side is increasing, it follows 


supE[r 2 £ s 4p ] < C [ su P E[r 2 £ 4p ]ds + CAt 4p , 

s<t Jo u<s 

from where we conclude the result thanks to Gronwall’s Inequality. 

Now we present the proof of the technical lemmas. 

Proof of Lemma 3.3. By the integration by parts formula. 


E[r 2 £ 4p ] = 4pE 


T 2 ^- 1 {sgn(Z 8 )&(*,m) ^ KX S )} ds 


+ 2p(4p — 1)E 

nt 


T 2 £ 4p ~ 2 E 2 ds 


■ 4pE 


U 0 


U o 


-E 


2fi s T 2 s £t P ds 


L^O 


□ 
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Thanks to Lemma 2.4 and the control in the moments of the exact process in Lemma 1.2 we have 

' T 2 s £t p - 1 dL° s (X) 


E 


< E 


| X 4p ~ 4 \dL° s {X) 


La. 

VI 

sup A s P 2 

V 

Lo<s<T 


E [L° (A') 2 ] < CAt 4p . 


On the other hand, with sgn(x) = 1 — 21r x<0 }, calling AIL, = 1L S — ILL,), we get for all 0 < s < t 

_^, 1 2 r _-i 1 2 


E^<2 


ctA“ - ctA, 


+ 2 


~^?OL -rra n~^2(X— 1 . ___ 

crX s - aX v{s) - OLO X v{s) AIL, 


+ Ra ^{Z s <0}’ 


where we put aside all the terms multiplied by 1 

Rf :=4 [< (3) + av*X%jAW a 
x { [aX“ (s) + c* CT 2 A X ) 1 AW S 


{Zs< 0 } 


in 


-a. rQ: 


aX:-aX; {s) -aa 2 X 2 v a is) 1 AW s 


aXf - aX 

So, from the previous computations, the Lipschitz property of b, and Young’s Inequality, we conclude 

E[r 2 t £t p ] < 4p f E [r 2 e £t p ~ 1 (&(*„«) - b(X s ))) ds 


}■ 


/0 

+4p\\b'\ 

Jo 

+Ap(Ap — 1)E 


E 


p2 c4p 
1 v(s) C rj(s) 


ds + ApWb'Woo / E[r% 4 P]ds 


^£ 4p ~ 2 ^ 


ctA“ - ctA 0 ds 


+ (4p — 2)(4p — 1)E 


r 2 £ 4p ds 


+2(4p-l) f E[D,(A') 4p ]ds 
Jo 


-E 


2f3 s T 2 s £ 4p ds 


UO 


E 


R^{Z b < 0} 


ds + CAt 4p . 


where R s = Ap(Ap — l)£ 4p 2 i?f + 8 p£ 4p 1 6(A I? (,)), and from Lemma 2.2 we have 


E 


R^{Z S < 0} 


< CAt 4p . 


Since Af < A max (a) and Remark 3.2-(*), we can apply Lemma 2.5 so E[D,(A) 4p ] < CAt 4p . Introducing these 
estimations in the previous computations, we have 


E[T 2 t £f p ] <4 p f* E [r 2 £ 4p_1 (b(X v{s) ) - &(A,))] ds 
Jo 


■4p||6'j 


E 


p 2 c4p 

1 v{s) C r](s) 


ds 


+ (^ll&'lloo + (4 p - 2)(4 p - 1)) E 


T 2 £ 4p ds 


Uo 


- 4p(4p — 1)E 


r 


Uo 


l £ 4 P -2 ^ 


crA'f -<rX.) ds 


- E 


2p s T 2 s £f p ds 


+ CAt 4p . 
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(3.7) 


Now we use the particular form of the weight process. Since for all \ < a < 1, for all x > 0, y > 0, 


x a -y a \(x 1 - a + y i - a )<2a\x-y\, 


we have 


E 


4 a 2 a 2 

x 2(i- a ) ds ’ 


t ^2 c^P 
1 v( s ) C ri{s) 


ds + CAt ip . 


; f (aXf - aX a s ) 2 ds < E f T 2 £j p - 

J o v ' Jo 

and then, from the definition of /3 in (3.1), we conclude 

E[r^ 4p ] < 4p f E [r% 4p - 1 (b(X v{s) ) - 6(X S ))] ds + 4p||& , || 00 f E 
Jo Jo 

from where 

E[T 2 t £t P ]<Ap /' 4 E[r^- 1 (6(X, (s) )-6(X s ))]ds + 4p||&'|| 00 f sup E[r 2 ^]ds + CAi 4 ?. 

Jo Jo 0<u<s 

Proof of Lemma 3.4. By integration by parts 

E [rfo 4p - 1 (&(* Ij(8) ) - 6(X S ))] = —E / 2r. u £: 4 P- 1 (6(X, )(s) ) - b{X u ))p u Y u du 

Jrt(s) 


□ 


E f r 2 d(C- 1 (6(X r?(s} )-6(X„))). 

Jri(s) 


(3.8) 


'V(s) 

Applying Holder’s Inequality to the first term in the right-hand side we have 


E [ rl^MX^) - b(X u )\/3 u du 

Jri(s) 


< / (E [r 2 ^]) 1 * (E [|6(X, W ) - b(X u )\ 4p pt P ]) 4p du. 


Recalling the Remark 3.2 -(ii), we have that E[/3® p ] is finite, so applying Lemma 2.1, 

E - b(X u )\ 4p /3t P ] < v /e[|KAV s) )-6(X u )|8p]E[/ 3® p ] < CAt 2p . 


Then, 


E f Tlf^HX^) - b(X u )]/3 u du <c(supE[r 2 5 4p ]) ^ At 3 / 2 . 

J r)(s) \u<s J 


Applying the Ito’s Formula to the second term in the right-hand side of (3.8), and taking expectation we get 

e f r 2 l d(£t p ~ 1 [b(x ri{s) )-b(x u ))) 


lri(s) 


= (b'(X u )b(X u ) + ^-b"(X u )X 2 ^ du 

+ (4 p- 1)E / r 2 f 4p_2 (6(X^ (s) ) - b(X u )) {sgn (Z u )b(X v{s) ) - b(X u )} du 

Jr){s) 

+ ^(4 P - 1)(4 p - 2)E f r 2 u £t p - 3 (b(Xv(s)) ~ b(X u ))Y 2 u du 
A Jri(s) 

-(4p-l)a 2 E f T 2 X p -' 2 b'(X u )X^ u du 

Jri(s) 

+ fc^E [ S T 2 u £t p - 2 {b{X^ s ))-b{X u ))dL 0 u {X) 

^ Jri(s) 


(3.9) 


’V{s) 

= '■ II + 12 + I 3 + 1*1 + h ■ 
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By the finiteness of the moment of X, the linear growth of b, and the polynomial growth of b ", applying Holder’s 
inequality, we have 


\h\<c (Efr^]) 1 4 * E 

Jv(s) 


\b'{X u )b{X u ) + -b"(X u )X 2 u ^ p 


du 


< C fsupE[r^P]) JP Af. 

\u<s J 

For the bound of / 2 > since b is Lipschitz, and sgn(:r) = 1 — 2tf x<0 \, we have 
|/ 2 | <C f E [T 2 X p - 2 \X v{s) - XnWX^ - X u \] du 

Jri(s) 

+ c f E[rl£t p - 1 \X v{s) -X u \\du+ f E\RWt { z u<0} \du 

Jri(s) "'fjO) 

<C fsupE[r2£4p] N ) _p Af 2 + C ("supE[r 2 £^ p ] N ) 4p At 3/2 + CAt 4p 

\u<s J \u<s J 

Where again, all the terms multiplied by l^ z <(| , are putted in the rest n„ , and the expectation of the product is 
bounded with Lemma 2.2. 

In a similar way for the bound of / 3 , decomposing E„ with sgn(a;) = 1 — 21{ x<0 }, 

|/ 3 1 <c f E [r 2 ^p-3| Xr)(s) _ Xu \D s (X) 2 ] du 

J ri(s) 


lri(s) 

+ C [ E 

Jri(s) 


T 2 u £t p - 3 \X V{S) - X u \ (aX u - aXZf\ du+ / E\R^\^ <Q} \du. 

J dn(s) 


For the first term in the right-hand side we have 

E [r 2 u £^ 3 \x v{s) -X U \D S (X) 2 ] < (E[r 2 C]) 1_ * (E[|X„ (8 ) -x u \ ip })^ (E[D S (X) 4 ^ 

1_3_ 

<(supE[r 2 C]) 4P At 5 / 2 , 

\u<s J 

due to the bound for the increments of the exact process and Lemma 2.5. For the second term, applying (3.7), and 
noting that from Remark 3.2 -(ii), the exact process has negative moments up to of order 12 p + 2 

E {T 2 u £* p - 3 \X vis) - X u \(aX u - aX“) 2 ] 

< CE\rl£?- l \X v(a) - X u \—p2—^] 

-X-u 

< C (E[r 2 f 4p ]) 1_ ^ {E[\x vis) - X u \ sp )) 1/8p (e 

1—— 


1 


X.. 


16(1 — a)p 


1/8 V 


< (su P E[r 2 ^]) 4p At 1 ' 2 . 

\u<s J 

We control the third term in the right-hand side in the bound for |/ 3 | using again Lemma 2.2, so 

|/ 3 | < C fsu P E[r 2 ^]) 4P At 7 ' 2 + fsu P E[r 2 ^]) 4P At 3 ' 2 + CAt ip . 

\U<S J J 
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Now we bound | J 4 1. 


\h\ < C f E\T 2 u E%>- 2 D u (X)b\X u )XZ\du 

Jri(s) 

+ c f nvl£t v -\aX u -aXZ)b'{X u )XZ\du+ f E\R^l { ^ <0} \du. 

Jri(s) d-n(s) 

We control the first term in the right-hand side using Holder’s inequality. Lemma 2.5 and the control in the moments 
of the exact process for all 0 < u < s 

E\r 2 u £t P ~ 2 D u (X)b'(X u )XZ\ < (E[r^])i-^(E[^(X)^])^(E[y(X u ) 4 ^“])4 

<C (sup[Tl£^]\ At. 

\u<s J 

For the second term in the right-hand side of the bound for |Z 4 1, we use one more time (3.7), and the existence of 
negative moments of the exact process X, and then 

E\T 2 u £t P ~ 2 [vX u - <] b'(X u )X :| < cnrlE^-^b'iXum < C (supEfT^]) ^ ■ 

Xu \u<s J 

To control the third term in the right-hand side of the bound for | Z 4 1 we use Lemma 2.2 just as before. So 

\h\ < c (supEpr 2 ^]) 2p At 2 + C (supEfT 2 ^]) 4p At + CAt 4p . 

\u<.s J \u<s / 


Finally, 


141 = f r2 u X u P ~ 2 \&(*„«) - b(X u )\dL° u (X) 


< CE 


’V{s) 

sup [1 + X^- 1 ] L° t (X) 

U<S 


< C \ E 


sup 

n<s L 


1 + xtr 1 


E[L°(A') 2 ] < CAt 4p , 


the last inequality comes from Lemmas 1.2 and 2.4. 
Putting all the last calculations in (3.8) we obtain 


1 T 


- b(X s ))}\ < C supE[r 2 ^] At 3 / 2 + c supEfT 2 ^] At 

\U<S J \U<S J 

1 _ 1 _ ^_ 3 _ 

+ c fsupEir 2 ^]) 2p Af 2 + c fsupE[r 2 C]) 4p At 7 ' 2 

\u<S J \U<S J 


CAt 4p . 


Applying Young’s Inequality in all terms in the right, we get the desired inequality (3.6). □ 

Remark 3.5. Proceeding as in the Proof of Lemma 3.4, ifp = 1 orp > 3/2, is not difficult to prove 

|E [£l p -\b{X v{s) ) - b(X s ))] | < CsupE^] + CAt 2p . (3.10) 

U<s 
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Notice that (3.10) is similar to (3.6) with the process F = 1. The restriction on p comes from the following observation. 
Applying the ltd’s Lemma to the function ( x,y ) K > x 2p ~ 1 y for p > 1, with the couple of processes — 

b(X.)), in the analogous of the identity (3.9) with T = 1, it will appear a term of the form 

eT £l p -*[b(X^ {s) )-b(X u )}Y? u du. 

Jn(s) 

The restriction p > 3/2 avoids the situation where 2p — 3 is negative. 

3.2 Proof of Theorem 1.6 

We start by controlling the L 2p -error. First, assume p = 1 or p> 3/2. By the Ito’s formula we have 


E[£ 2p )=2pE f* £ 2p - 1 {sgn(Z s )b(X v(s) )~b(X s )}ds + 2pE f T^dL^X) + p(2p - 1)E f £ 2p ~ 2 E 2 s ds. 
Jo Jo Jo 

As we have seen before, E f* £ 2p ~ 1 dL°(X) < CAt 2p , and sgn(x) = 1 — 21{, E< o}, so 

E[£ 2p ] < 2pE f £?-' [b(X v{s) ) - b(X v(s} )] ds 
Jo 

+ 2 P E f T 2 /- 1 [&(X„ (a) ) - b(X s )] ds 
Jo 

+ 8p(2p - 1)E [ £ 2p ~ 2 Lx: ~ 

Jo 

+ 8p(2p — 1)E^ £ 2 s p ~ 2 D s (X) 2 ds + E R s t { z s<0} ds + CAt 2p , 


(3.11) 


ds 

rt 


where R s = Ap(2p - l)£ 2p 2 + ao 2 X 2 r ^ 1 (W s - lF/(s)) + 1 ^(^( s ))- If we use the Lipschitz 

property of b, and Young’s inequality in the first term in the right of (3.11), Lemma 2.5 in the fourth one, and 
Lemma 2.2 in the fifth one, we have 


E[£ f 2p ] < <7 / su P E[£ 2p ]ds + 2 p / E[£ 2p-1 (&(X )?(s) ) - b(X s ))]ds 


/0 u<s 


+ 8p(2p — 1) f E[£ 2 s p ~ 2 (aX a a ~ oX?) 2 ]ds + CAt 2p . 
Jo 


(3.12) 


And, according to Remark 3.5, we have 


|E- 6(X a ))]| < CsupE^P] + CAt 2p . 

IL<S 


On the other hand, using again (3.7), we have 

e [£ 2 s p ~ 2 (oX a s - ax/) 2 ] < ce\£: p X f 2{1 ~ a) ] = ce.\t B e^ 


and applying Cauchy-Schwartz inequality. 


E[r s ^- 


i 


2(1 —a) 




r/ 1 ] <(E[rX 4p ])= 



1 

2 


The first term in the right-hand side is the weight error controlled by Proposition 3.1. To control the second one, let 
us recall Remark 3.2. For a >b the exact process and the weight process T have negative moments of any order. 
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therefore the second term in the last inequality is bounded by a constant. On the other hand, for a = ^ we need a 
finer analysis. From the second point in Remark 3.2, the 12 p + 2-th negative moment of the exact process is finite, and 
since 

16p(4p - 1) 2 6?J + 1 < (12 p + 2) 2 , 
bp 

according with the third point of Remark 3.2, the 2(6 p + l)/6p-th negative moment of the weight process T is also 
finite. Therefore, when a = 


E 


1 

-^-4(1 -a) 




< 


E 


1 

j£-12p+2 


( 

o6p+l “ 

IE 

P 6p 

1 S 


6p 

6p+l 


<C, 


and then in any case 


E 



1 

^2(l-a) 



< CAt 2p . 


Introducing all the last computations in (3.12) we get 


E[£ 2p ] < C 


[ (supE [f 2p ]) ds + CAt 2p . 
Jo \u<s J 


Since the right-hand side is increasing, thanks to Gronwall’s Inequality we have, for p = 1 or p > 3/2, 

sup (e[£ 4 2p ]') 2p < CAt. 
o <t<T k ' 


We extend to p £ (1,3/2), thanks to Jensen’s inequality 


(E < (E[£ t S D f < CAt, 

and we conclude that (E[£ t 2p ])5p < CAt for all p > 1 satisfying Remark 3.2- (i). 
Now we control the //'-error. For o- = / and p > 2, denoting p' = p > 1, 

i 

sup (E[£’f])p= sup (E[£t P ]) 2p <CAt. 

o <t<T o <t<T k > 


Hypothesis 1 ,5-(z) gives 


6 ( 0 ) > 


3(2 p + l)cr 2 
2 


3(4 pf + l)a 2 

2 


Since (3.4) in Remark 3.2 is satisfied, we can control the L 2p (O )-norm of the error and then 


sup (E[£ p ]) p = sup (e[£ 2p ]V P < CAt. 

0 <t<T 0 <t<T k ) 


If p € [1,2), Hypothesis 1.5-(i) is 6(0) > 15^-, which is enough to bound the L 2 (fl)-norm of the error, and then 
from Jensen’s inequality 

sup (E[£ p ])p < sup (E[£ 2 \y <CAt. 

0 <t<T 0 <t<T 

The case a > 1/2 is easier. Since the Hypothesis in the parameters for this case does not depend on p, we can conclude 
for any p > 1 from Jensen’s inequality and the control for the T 2p (fl)-norm of the error. 
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Remark 3.6. Let us mention an example of extension of our convergence result, based on simple transformation 
method: consider the 3/2-model, namely the solution of 

r t = r 0 + [ cir s (c 2 - r s )ds + f c 3 rl ,2 d.W s . 

Jo Jo 

Applying the ltd's Formula to v t = f(r t ), with f(x) = rr , we have 

vt = v 0 + / (ci + C3 - cic 2 v s ) ds + / c 3 v^ 2 dB 8i 
Jo Jo 

where B s = — W s is a Brownian motion. We can approximate v with the SMS v, and then define f t : - 1 /v t . Then we 
can deduce the strong convergence with rate one of ft to rtfrom our previous results. 

Transformation methods can be used in a more exhaustive manner, in the context of CEV-like SDEs and we refer 
to [6] for approximation results and examples, using this approach. 


3.3 Strong Convergence of the Projected Milstein Scheme 

The Projected Milstein Scheme (PMS) is defined by 


2 \r2ot — l 


Xt k = 


(X_ x + b{X th _ x )M + aXZ.AWt, - W tk _0 + aa ^ [{W th - W th _0 2 - At] ) + , A 0 = cco 


where for all x £ R, ( x) + = max(0, x). The continuous time version of the (PMS) is given by 

T 2 X 2a ~ l 


X t — f(t) + b(X v (t))At + aX^(W t - W v{t) ) + 


■nit) 


[(Wt - W v{t) ) 2 - At] 


(3.13) 


Notice that for all t € [0, T\, 0 < X t < X t , then the positive moments of the PMS are bounded (see Lemma 1.4). 

To obtain a strong convergence rate for the PMS, we first show that the PMS and the SMS coincide with a large 
probability. 

Lemma 3.7. Let us consider the stopping time r = inf{s > 0 : X s 7 ^ X s }. Assume that 6(0) > 2a(l — a) 2 a 2 and 
At < 1/(2 K(a)). Then for any p > 1, 

P(r < T) < CAt p . 


Proof. Notice that r is almost surely strictly positive because both schemes start from the same deterministic initial 
condition xq. On the other hand 

JV-l 

P (r < T) = ^ P ( r e (*fc> 4+i]), 

2—0 


and according to Lemma 2.2 


P(t£ [ti,t i+ 1)) = P 


inf Z s < 0,X th 

tk |_i 



< 


inf Z s -- 

tk 


< 0 



So, 

p(T<r)<L exp (^_I). 

Since for any p > 1, there exists a constant C p such that exp(—j/At)/At < C p At p , we have 


P(r < T) < CAt p . 


□ 
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Corollary 3.8. Assume Hypotheses 1.1 and 1.5. Consider a maximum step size A max (a) defined in (1.10). Let 
(Xt, 0 < t < T) he the process defined on (1.1) and ( X t , 0 < t < T) the Projecter Milstein scheme given in (3.13). 
Then for any p > 1 that allows Hypotheses 1.5, there exists a constant C depending on p, T, 6(0), a, a, K, and Xq, 
but not on At, such that for all At < A max (o;), 

sup (E[\X t - X t \P]) 5 < CAt. (3.14) 

0 <t<T v ' 

Proof. Notice that for all t £ [0, T\, with r = inf{s > 0 : X s ^ X s }, 

E[|A t - X t \ p } = E[|X t - A t | p l {T < T} ] +E[|A t - X t \ p t {T>T} ] 

< ^E[\X t ~X t \ip]P(T<T) + E[|A t - X t \P] < CAt p 

where the last inequality comes from Lemma 3.7 and Theorem 1.6. □ 


4 Numerical Experiments and Conclusion 

We start this section with the analysis of two numerical experiments. The first one aims to study empirically the strong 
rate of convergence of the SMS in comparison with other schemes proposed the literature. The second one aims to 
study the impact of including the SMS in a Multilevel Monte Carlo application. 


4.1 Empirical study of the strong rate of converge 

In this experiment we compute the error of the schemes as a function of the step size A/ for different values of the 
parameters a and <j. 

For a > -5 we compare the SMS with the Symmetrized Euler Scheme (SES) introduced in [3], and with the 
Balanced Milstein Scheme (BMS) presented in [13]. Whereas for a = |, in addition to the aforementioned schemes, 
we will also compare SMS with the Modified Euler Scheme (MES) proposed in [ 6 ], and with the Alfonsi Implicit 
Scheme (AIS) proposed in [1]. 

Let us first, shortly review those different schemes. 


Alfonsi Implicit Scheme (AIS). Proposed in [1], the AIS can be applied to equation (1.1) when the drift is a linear 
function. A priori, the AIS can be applied for a £ [i, 1), but is relevant to observe that only when a = \, the AIS is 
in fact an explicit scheme (also know as drift-implicit square-root Euler approximations ) whereas in any other case is 
not. This implies that in order to compute the AIS for a > at each time step it is necessary to solve numerically 
a non linear equation. This extra step in the implementation of the scheme brings questions about the impact of the 
error of this subroutine in the error of the scheme, and about the computing performance of the scheme. Since this 
questions are beyond the scope of the present work, we include the AIS in the comparison only in the Cox-Ingersol- 
Ross (CIR) case (linear drift and a = ^ j. In this context, the AIS can be use only if cr 2 > 46(0), for other values of 
the parameters the AIS is not defined. In terms of convergence, when a = according to Theorem 2 in Alfonsi [2], 
the AIS converges in the L p (fl)-norm, forp > 1, at rate At when (1 V 3p/4)a 2 < 6(0). When a > ^, the AIS (see 
Section 3 of [2]) converge as soon as 6(0) > 0, at rate At to the exact solution. 


Balanced Milstein Scheme (BMS). The BMS was introduced by Kahl and Schurz in [13], and although its con¬ 
vergence it is not proven for Equation (1.1) (see Remark 5.12 in [13]), numerical experiments shows a competitive 
behavior (see [12]). Also, the BMS can be easily implemented for a £ [|, 1), so we decide to include it in our 
numerical comparison. 
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Modified Euler Scheme (MES). Introduced in [6], the MES can be applied to the Equation (1.1) for a £ [|, 1) 
when the drift has the form b(x) = Hi( x ) — H2(x)x for /ii and H2 suitables functions. The rate of convergence in the 
L 1 (fl)-norm of the MES depends on the parameters. For a = ^ the rate is 1 if a 2 is big enough compared with 6(0), 
and it is p < 1 in other case. When a > the MES converges at rate 1 as soon as b{ 0) > 0 (see Proposition 4.1 in 
[6]). When a > the implementation of the MES requires some extra tuning which is not explicitly given in [6] (see 
Remark 5.1), so we implement the MES only for a = 

Symmetryzed Euler Scheme (SES). The SES, introduced in [3], is an explicit scheme which can be apply to the 
equation (1.1) for a £ [ 5 , 1 ) and any Lipschitz drift function b. It has the weakest hypothesis over b of all the schemes 
discussed in this paper. If a = |, according to Theorem 2.2 in [3], the rate of convergence of the SES is \/At under 
suitable conditions for 6(0), a 2 and I\. When a> \ the SES converge at rate \JAt as soon as 6(0) > 0 (see Theorem 
2.2 in [3]). 

We summarize the theoretical analysis of the schemes above in Table 1 for a = ^, and in Table 2 for a > 


Scheme 

Norm 

Drift 

Convergence’s Condition 

Theoretical 

rate 

SMS 

LP,p > 1 

6 Lipschitz, b £ C 2 
b " with polynomial growth 

6(0) > 3 (2[p V 2] + 1) y 

1 

AIS [2] 

LP,p£[l,*W i) 

b(x) = a — bx 

6(0) > (1 V | p)a 2 

1 

BMS 




undetermined 

MES [6] 

L 1 

b(x) = Hi{x) - p, 2 {x)x 

Hi £ C 2 n , Hi > 0 

Hi < 0, H 2 > 0 

6(0) > 

1 

6(0) > ^l 2 

1 

2 

6(0) > fj 2 

7 1 1 

(6’ 2 2b(0)+a 2 J 

SES [3] 

LP,p> 1 

6 Lipschitz 

m> [V^(§vi) + iJ 
lC(q) = K(Wq - 1 ) 

V4o- 2 (8p- l) 2 

1 

2 


Table 1: Summary of the condition over the parameters for the convergence of the different schemes for a = ^. 


Simulation setup 


In our simulations we consider a time horizon T = 1, and xq = 1. In order to include as many schemes as possible 
we consider for all simulations a linear drift 

b(x) = 10 — 10x. 

To measure the error of each scheme, we estimate its L 1 (fl)-norm for which a theoretical rate is proposed for all the 
selected schemes. 

Let E|£™ s |, E|£™ s |, E|£^ ES |, E|£“ ES |, and E|£" s | be the L [ (O)-norm of the error for the SMS, BMS, SES, MES 
and AIS respectively. To estimate these quantities, we consider as a reference solution the AIS approximation for 
At = A max (a)/2 12 when a = and the SMS for At = A max (a)/2 12 when a > |. Then for each 


At £ 


( A max (a) 
\ 2 ™ 



we estimate E|£)V| by computing 5 x 10 4 trajectories of the corresponding scheme, and comparing them with the 
reference solution. The results of these simulations are reported in Figures 1 (a = and 2 (a > 1). The graphs 
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Scheme 

Norm 

Drift 

Convergence’s Condition 

Theoretical 

rate 

SMS 

LP,p> 1 

6 Lipschitz, 6 £ C 2 
b" with polynomial growth 

6(0) > 2a(l — a) 2 cr 2 

1 

AIS 

i p ,pe[i,f) 

b(x) = a — bx 

6(0) > 0 

1 

BMS 




undetermined 

MES 

L 1 

b{x) = pi{x) - p. 2 {x)x 

Pi £ c 2 nc£, pi > 0 
p[ < 0, p' 2 > o 

6(0) > 0 

1 

SES 

L p ,p > 1 

6 Lipschitz 

6(0) > 0 

1 

2 


Table 2: Summary of the condition over the parameters for the convergence of the different schemes when a > \. 


plot the LogE|£y'| in terms of the Log A t, and we have added the plot of the identity map to serve as reference for 
rate of order 1. The schemes with a slope smaller than the slope of the reference line have an order of convergence 
smaller than one. To obtain a more quantitative comparison of the schemes, we also perform a regression analysis on 
the model 

log(E|fj' I) = plog(Af) + C. 

Notice that p. the estimated value for p, corresponds to the empirical rate of convergence of the different schemes. We 
present the result of this regression analysis in Tables 3 and 4. 

Empirical results for a = Figure 1 and Table 3 present the result for the CIR case. From Table 1, we observe 
that we can distinguish five cases for the parameters. 

The first case (a 2 = 1) is such that 6(0) > 6cr 2 : the SMS, the AIS, and the MES have a theoretical rate of 
convergence equal to At, whereas the SES has a theoretical rate of convergence equal to V At. In Figure la, we 
observe that the graphs of the SMS, the AIS, the BMS, and the MES seem parallel to the reference line, which is 
expected, while the SES has a smaller slope. This is also confirms in the first line of the Table 3, where we observe 
that the empirical rates of convergence are close to the theoretical ones. Notice that the BMS has a competitive 
empirical rate of convergence, although the theoretical one is not known. 

The second case (u 2 = 4) is such that 6(0) £ (5^-, 6er 2 ), now only the AIS and the MES have a theoretical rate 
of convergence equal to At. However, how we can see in Figure lb the SMS still shows a linear behavior in this 
case. Recall that the condition over the parameters is a sufficient condition and we believe that could be improved. 
Notice that also the BMS shows a linear behavior. In the second line of Table 3, we can observe that empirical rates 
of convergence are close to one for all the scheme but the SES. 

In Figure lc, we illustrate the third case (ex 2 = 6.25) and then 6(0) £ (3cr 2 /2, 5 ct 2 /2). In this case, only the AIS 
has a theoretical rate of convergence equal to one. For the MES is y/At, but in the graphics we still observe a linear 
behavior for the MES, and also for the SMS and the BMS. This is confirm in the third line of Table 3. 

The fourth case is (a 2 = 9) and 6(0) £ (cr 2 ,3 <t 2 /2), which we display in Figure Id. For this values of the 

parameters the theoretical rate of convergence is known only for the AIS and the MES. Nevertheless, we observe in 
the graphs and in the fourth line of Table 3 that all the schemes seems to reach their optimal convergence rates. 

Finally, the fifth case is (cr 2 = 36) and then 6(0) < a 2 . In this case all the schemes have a sublinear behavior 
as we can see in Figure Id and the fifth line of Table 3. This case illustrate the necessity of some condition over the 
parameters of the model to obtain the optimal rate of convergence for the SMS. 

Empirical results for a > ^. In Figure 2 and Table 4 we present the results of the simulations for a = 0.6, and 

a = 0.7. 
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Observed L l (Vl] 

convergence rate 6 (and its R 2 

value) 



a 2 

SMS 

AIS 

BMS 

MES 

SES 

P 

C R 2 ) 

P 

(R 2 ) 

P 

(R 2 ) 

P 

(R 2 ) 

P 

C R 2 ) 

1 

0.9956 

(99.9%) 

1.0060 

(99.9%) 

1.0055 

(99.9%) 

0.9955 

(99.9%) 

0.5941 

(99.3%) 

4 

0.9969 

(99.9%) 

1.0054 

(99.9%) 

1.0037 

(99.9%) 

0.9961 

(99.9%) 

0.5344 

(99.8%) 

6.25 

0.9976 

(99.9%) 

1.0043 

(99.9%) 

1.0002 

(99.9%) 

0.9941 

(99.9%) 

0.5237 

(99.9%) 

9 

0.9984 

(99.9%) 

1.0015 

(99.9%) 

0.9859 

(99.9%) 

0.7891 

(99.9%) 

0.5164 

(99.9%) 

36 

0.6410 

(99.7%) 

0.6282 

(99.8%) 

0.4538 

(99.3%) 

0.3575 

(99.4%) 

0.4718 

(99.9%) 


Table 3: Empirical rate of convergence p for the L 1 (fi)-error of the schemes when a = \ for different values of a 2 . 


In these cases, it can be observed in numerical experiments that the MES needs smaller At to achieve its theoretical 
order one convergence rate, unless one tunes the projection operator in the manner of Remark 5.1 in [6]. Since this 
tuning is not explicitly given we do not include the MES in these simulations. 


Parameters 

Observed L x ( ft) convergence rate p (and its R 2 

value) 



SMS 

BMS 

SES 

a 

a 2 

P 

(R 2 ) 

P 

(R 2 ) 

P 

C R 2 ) 


49 

0.9819 

(99.9%) 

0.7296 

(99.1%) 

0.5273 

(99.8%) 

0.6 

53.29 

0.9766 

(99.9%) 

0.7788 

(99.3%) 

0.5133 

(99.9%) 


144 

0.6609 

(98.9%) 

0.4336 

(97.3%) 

0.5074 

(99.9%) 


64 

1.004 

(99.9%) 

0.9022 

(99.7%) 

0.5242 

(99.8%) 

0.7 

81 

0.9991 

(99.9%) 

0.8813 

(99.7%) 

0.5327 

(99.7%) 


225 

0.9146 

(99.7%) 

0.6497 

(97.6%) 

0.6410 

(99.2%) 


Table 4: Empirical rate of convergence p for the 1? (O)-error, when o: > | for different values of a and a 2 . 

We have observed in the numerical experiments three cases for the parameters. The first one is when 6(0) > 
2o( 1 — a) 2 cr 2 (a 2 = 49 and a 2 = 64). In this case. Theorem 1.6 holds and we observe the order one convergence 
(see Figures 2a, 2b, and first and fourth row in Table 4). The second case is when the parameters do not satisfy 
6(0) > 2a(l — a) 2 a 2 ( a 2 = 53.29 and a 2 = 81), and then we can not apply Theorem 1.6, but in the numerical 
simulations we still observe the order one convergence (see Figures 2c, 2d, and second and fifth row in Table 4). 
Finally the third case, is when a 6(0), and then we do not observe a linear convergence anymore (see Figures 2e, 
2f, and third and six row in Table 4). Notice that in the three cases the SMS performs better than the BMS, specially 
when <t 2 grows, (see Table 4). 
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The second and third case show us that some restriction has to be impose on the parameters to observe the con¬ 
vergence of order one. But our restriction, although sufficient, it seems to be too strong, specially for a close to 
one. 


4.2 Application of the SMS in Multilevel Monte Carlo 


We continue this section by testing the SMS in the context of a multilevel Monte Carlo application widely used 
nowadays in computational finance (see e.g. [11] and references therein). Multilevel Monte Carlo is an efficient 
technique introduced by Giles [9] to decrease the computational complexity of an estimator combining Monte Carlo 
simulation and time discretisation scheme for a given threshold in the accuracy. For details we refer to [9, 10, 1 1], 

For this experiment, we consider the classical but non trivial test-case of the Zero Coupon Bound (ZCB) pricing 
of maturity T, 


B(0,T) = E 



under the hypothesis that the short interest rate dynamics ( r t ,t > 0) is modeled with a CIR process (a = A and 

b{x) = a — bx ): 

dr t = (a - br t )dt + a^ t dW t . 


In this context, the price of the ZCB admits a wellknown closed-form solution given by (see e.g [7, 14]) 


B(0,T) = A(T)e" B(T)ro , 


A{T) = 


2 Ae (6+A)T/2 

(A + b)(e XT — 1) + 2A 


B(T) 


2 (e AT - 1) 

(X + b){e XT - 1) + 2A' 


where r 0 is the initial value of the interest rate, and A = y/b 2 + 2a 2 . Let EB(At^) a discrete-time weak approxima¬ 
tion of B( 0, T) with the time step A t^y We consider the /.-level Monte Carlo estimator : 


-j Nq L Ni 

Yt = w w £ (b“>(a t m ) - b«>(a 


, . N t 
1 = 1 2=1 


For a targeted mean-square error e 2 on the computation of the quantity //((), T) 


E[(Y T -B(0,T)) 2 } = O(e 2 ), 


one can choose the following a priori parametrization of the MLMC method in order to minimize the computational 

log e -1 

time (complexity) (see [9, 10, 11]): we use the estimation L = -; from one level to the next one the time step is 

log 2 

divided by 2, A - 2 ,,^ l , ; the number of trajectories to simulate is estimated with Giles formula [9] 


Ni = ^y/ViAti (j2yJVi/At {l) 

\l =0 

with Vi = Y cir^B^ (A t^) — B ^ (A■ As an estimator for the bias variance, the strong rate of convergence 
of the discretisation scheme enters as a key ingredient in the TV; a priori estimation. A scheme with a reduced strong 
bias will then allow a smaller TV;. We apply the MLMC computation for the SMS, the PMS, the AIS and the BMS. 

We summarize the results of the performance comparison between the four schemes in Table 5. The computation 
have been run using a initial interest rate r 0 = 1, the maturity of the bond T = 1, the drift parameters a = h = 10, 
the volatility cr = 1. For the MLMC simulation, we fix a minimum number of trajectories equal to 500, and a minimal 
number of levels equal to 6. 

In Table 5, we give the measures of the CPU time for a set of three decreasing targeted errors, as long as the 
effective measured error and the total number of simulated trajectories. As expected, the required threshold error has 
been reached by the MLMC strategy. As also expected (see Giles [10]), Milstein schemes perform better than their 
Euler versions. Finally, as e decreases, the SMS clearly performs better than his PMS version. 
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log {Error) log {Error) 




(a) Parameters in case 1: 6(0) > 6cr 2 . (b) Parameters in case 2: 6(0) G (5 ct 2 /2, 6<t 2 ). 




(c) Parameters in case 3: 6(0) G (3<t 2 /2, 5ct 2 /2). (d) Parameters in case 4: 6(0) G (a 2 , 3<t 2 /2). 



(e) Parameters in case 5: 6(0) < a 2 . 

Figure 1: Step size At versus the estimated I? (O)-strong error for the CIR Process (Log-Log scale). The identity map 
serves as a reference line of rate one. 
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log {Error) log {Error) log (Error) 




(a) Parameters in case 1: a = 0.60. (b) Parameters in case 1: a = 0.70. 




(c) Parameters in case 2: a = 0.60. (d) Parameters in case 2: a = 0.70. 




(e) Parameters in case 3: a = 0.60 (f) Parameters in case 3: a = 0.70 

Figure 2: Step size At versus the estimated L ] (O)-error for a > \ and different values for tr 2 (Log-Log scale). The 
identity map serves as a reference line of rate one. . 
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e = 1.0e-03 

(L = 9, A t (L) = 1/2 10 ) 

SMS 

PMS 

AIS 

BMS 

CPU time 

0.2304 

0.2657 

0.264 

0.274 

(JVo + • • • + N l ) 

(observed error) 

(792 651) 
(1.970e-05) 

(950 838) 
(3.347e-04) 

(990 769) 
(3.132e-04) 

(992 432) 
(3.292e-04) 

e = 1.0e-04 

(L = 13, A t (L) = 1/2 14 ) 

SMS 

PMS 

AIS 

BMS 

CPU time 

16.871 

20.843 

17.311 

16.95 

(No H-+ N l ) 

(observed error) 

(56 229 224) 

(4.870e-05) 

(70 876 600) 
(1.091e-04) 

(73 824 621) 

(9.538e-06) 

(73 668 115) 

(2.203e-06) 

e = 1.0e-05 

(L = 16, A t (L) = 1/2 17 ) 

SMS 

PMS 

AIS 

BMS 

CPU time 

1589.6 

1910.7 

1576.4 

1540.2 

(No H-+ N l ) 

(observed error) 

(5 531 879 264) 

(4.752e-06) 

(6 913 546 698) 

(5.889e-06) 

(7 368 734 119) 

(3.912e-06) 

(7 333 474 098) 

(5.653e-07) 


Table 5: CPU time to achieve the target error for the different schemes. The observed error is | Yt — B( 0, T)\. 


4.3 Conclusion 


In this paper we have recovered the classical rate of convergence of the Milstein scheme in a context of non smooth 
diffusion coefficient, although we have to impose some restrictions over the parameters of the SDE (1.1) to ensure 
the theoretical order one of convergence. Typically, if the quotient 6(0) /a 2 is big enough we will observe the optimal 
convergence rate. 


In the numerical simulations we have observed that, despite the fact it is necessary to impose some restriction over 
the parameters of SDE (1.1) to obtain the order one convergence. Hypothesis 1.5 seems to be not optimal, specially for 
a = Also, through numerical simulations, we have observed that the use of SMS could improve the computation 
times in a Multilevel Monte Carlo framework, at least as well as the (CIR specialized) one-order schemes. 


Although our result seems more restrictive in term of hypotheses on the set of parameters, in particular if we 
compare SMS with Lamperti’s transformation-based schemes (see the recent works in [2] and [6], the SMS can be 
applied to a more general class of drifts functions and in various contexts. It is thus a useful complement of the existing 
literature. 
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5 Proofs for preliminary lemmas 

5.1 On the Local Error of the SMS 

Proof of Lemma 2.1. From the definition of X, and the algebraic inequality for positive real numbers (a \ + ... + 
a n ) p < n p (a p + ... + a p ) we have 

\*t ~ X v{t) \ 2p < 3 2p (b(X v(t) ) 2p (t - p(t)) 2p + a 2p xlw(W t - W r)(t) ) 2p 

+ a S?4 2 ( r )2p [w w vw) 2 - - vm 2p ). 

Thanks to the linear growth of b. Lemma 1.4 and the properties of the Brownian Motion it is quite easy to conclude 
the existence of a constant C such that 

E[|X t -X„ (t) | 2p ] < CAt p , 


from where the result follows. 

5.2 On the Probability of SMS being close to zero 

From b (7 {a) and K(a) defined in (1.9), let us recall the notation 

M a ) 


□ 


=(a) := 


K(a) 


introduced in Lemma 2.3. As 6 , 7 ( 0 :) > 0 under Hypothesis 1.5-(«), x(a) is bounded away from 0. In particular, 

,. , (6(0)-(x 2 /4) u , 6(0) 

hm x(a) = -—-, whereas hm x(a) = — - 5—. 

, 1 A a-¥ 1 A + <j A 2 

“^2 ' 

Proof of Lemma 2.3. Denoting A W s = (W s — W^tX), and As = s — r?(s), we have for all s € [0, T], 

—= OCT" —2a— 1 . __ _o a a ,, - -r— (, \ OLO~ —2a- a , . 

Z s = —^( s ) XW S +aX r) ^XW s + A^( s ) + |^6(X^( S )) —X^ ) As. 

From the Lipschitz property of 6 and the following bound for any x > 0 


2(l-o) 


we have 


< 4(1 - a) 2 + (2a - 1)[2(1 - a)]^i, 


> ^X^AW 2 + aX; (s) AW s + X vW + (b a (o) - K(a)X v(s) ) As. 


(5.1) 


(5.2) 


So, 


Z s < (1 - p)b a (a)As, X v ( s ) < px(a) 
.2 


< 


~ 2 ~-^"n(s) XW 2 + aX^^AWs + X v ^ + (pb a (a ) — K(a)X v ^) As < 0, X v ^ < px(a) 

From the independence of A W s with respect to if we denote by A" a standard Gaussian variable, we have 

P Z s < (1 - p)b a (a)As 1 X v ( s ) < px(a) 


< E 


x2a + <x^Asx a ff + x + [pb a (a) — A'(o)a;] As < 0 


V — ^77(s) J 
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Notice that in the right-hand side we have a quadratic polynomial of a standard Gaussian random variable. Let us 
compute its discriminant: 

A(x, a) = a 2 x 2a As — 2aa 2 x 2a ~ 1 As (x + (pb r7 (a) — K(a)x) As) 

= —( 2a — l)a 2 x 2a As — 2aa 2 x 2a ^ 1 As 2 ( pb a {a ) — K(a)x). 

Since b cr (a) > 0, we have A(x, a) < 0 for all a £ [|, 1), and x < px(a). So, for all As < At we have 

P\Z S <(1- p)b a {a)As, X V ( S ) < px(a)] = 0, 

taking As = At we conclude on the Lemma. □ 

Proof of Lemma 2.2. We have 


inf Z s < 0 ) =: 


tk<S<tk +1 


inf Z s < 0,X tk > x(a)\ +P [ inf Z s < 0, X tk < x(a) ) (5.3) 


tk<S<tf E + i 


We start with the second term in the right hand of the last inequality. By continuity of the path of Z and Lemma 2.3, 
we have 

P( inf Z s < 0, X tk < x(a) ) = Y' P (Z s < 0, X tk < x(a)) = 0. 
v ' seQn(t fc ,tfe+i] 

On the other hand, from 5.2 we have 

Z s > oX a v{s) AW s + (1 - K(a)As)X vW + b a {a)As. 


(5.4) 


Then 


inf Z s < 0, X tk > x(a 

t k <S<t k+ 1 


< P ( inf 


■Z r/(s) ( b a (cx ) K (cr) A 77 ( 5 )) As 


oX 


+ A1L S < 0, X tk > x(a) 


vM 


= E fj{X tk ) l{x tjl >g(a)}_ > 

where the last equality holds thanks to the Markov Property of the Brownian motion, for 

x x ~ a b a (a) — K(a)x 


iP{x)=p( inf -- + - [a ’^u + B u <o), 

v ' \o<u<At cr ax a > 


where (B t ) denotes a Brownian Motion independent of ( W t ). 

If (^,0 < t < T) is a Brownian motion with drift p, starting at yo- then for all y < yo , we have (see [4]): 


inf B£ < y 
o<s<t 3 ~ y 


= 1 erfc ( y^Zl + rff\ 

V2t y/2 J 


+ 2 ex P (~ 2 P-(yo - y )) erfc 


yo-y 


psjt 


Y 2 1 V 2 J ’ 


(5.5) 


where for zfM, erfc z = \/ 2 / 7 r exp (— u 2 /2) du. In our case y 0 = a: 1 a /a, p = (b„(a) — K(a)x) / ax a , and 
y = 0. Then 


i*x) = ierfc f 1(1 ~ A»Ai)i + MgjM \ 
2 V \/2Atcrx a J 

2[b a (a) — K(a)x]x 


- exp - 


erfc 


x — [b a (a ) — K(a)x]At 
\j2At<jx a 
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Since At < 1/2 K(a), for any x > x(a), the arguments in the erfc function in the last equality are both positives, and 
then recalling that for all z > 0 erfc(2) < exp(— z 2 ), we obtain 


Mx] < I exp f [(1 ~ K{a)At)x + b a (q)Atf 
nX) S 2 6 P ^ 2A ta 2 x 2 <* 

1 ( 2[b cr (a) — Jv(a)a:]aA 

+ 2 exp (.-^-J«P 

2 ' 


[x — [fecr(a) — K(a)x\At]- 
2Ata 2 x 2a 


< exp — 


[(1 — K(a)At)x + b a (a)At\ 


2Ata 2 x 2a 


So, for all x > x(a) 


ip{x) < exp ( — 


(1 — K(a)At) 2 x 2 ^ 1 “■* 
2a 2 At 


Then 


■( 


inf Z s < 0, X tk > x(a )) < E 


< exp ( — 

and finally, choosing 7 = a;(a) 2 ( 1_a ) /8a 2 , we get 

pf inf Z s < 0) < exp (— . 

\tk<s<tk +1 / V A t/ 


“ P { - 2^At - ) V-..>•<«» 

(1 — K(a)At) 2 x(a ) 2 ( 1 ~ a ) 

2a 2 At 


□ 


5.3 On the Local Time of the SMS at Zero 
The Stopping Times (0 Q , \ < a < 1) 

In what follows, we consider 

Q a = inf {s > 0 : X s < (1 — y/a)b a (a)At} . (5.6) 

Lemma 5.1. Assume b(0) > 2a(l — a) 2 a 2 , and At < \/{2K(a)) Aato/[(l — y/a)b a (a)}. Then there exists a positive 
constant 7 depending on a, b( 0), K and a but not on At such that 

ne„<T)<l t ^(X). (5.7) 

Proof. First, notice that the condition At < Xq /[{1 — -Ja)b a (a )] ensures that the stopping time 0 Q is almost surely 
strictly positive. 

To enlighten the notation along this proof, let us call l a (a) := (1 — y/a)b a (a), and C,k = inft fc < s <t fc+ i Z s . We 
split the proof in three steps. 

Step 1. Let us prove that for a suitable function xp : ffi. —► [0,1] and the set Ak = {Xt k > x(a)y/a} £ Tt k '. 


N -1 

P(Q a <T)< J2 K (^t k )t Ak ) 


k =0 


(5.8) 
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Indeed, 


N—l 


’(©« < T) < Y, p (Cfc < la(a)At, X tk > l„(a)At). 


k=0 


But, for each k = 0,..., N — 1 

P (Cfc < la (a) At, X tk > l a (a)At) = P (Cfc < l a (a)At, X tk > l a (a)At, X tk < x(a)^cx) 

+ P (Cfc < la(a)At, X tk > l a (a)At, X tk > x(a)y/a) . 

Since l a {a)At = (1 — yJa)b a {a)At < x(a)y/a, we have 

p(Cfc < iff (a) At, X tk > l a (a)At, X tk < x(a)y/aj 

< p (Cfc < iff (a) At, X tk < x(a)y/a) 

< ^2 P (Z s < iff(a)At, X tfc < *(a)\/a) = 0, 

seQn(t fc ,t fc+ i] 

thanks to Lemma 2.3. On the other hand, we have 

p(Cfc < iff(a)At, X tfc > l a (a)At, X tk > x(a)y/cxj 

= P (Cfc < iff(a)At, A tfe > x(a;)v / a) 




inf 


-=-=■1 — Q; 


^v( s ) (pa {a) K (a)A'.^( s )) As , l a (a)At , /— 

-l-AVL s <^=4—, A tfc > x(a)Va 


ctA 


?)(s 


crA 


7j(s) 




where the inequality comes from (5.2), and the last equality holds thanks to the Markov Property of the Brownian 
motion, for 


ip(x) =P( c 


x la b er (a)-K(a) x l a {a)At 

mt-1 - u + t> u < - 


, 0 <u<At <J <JX a 

where (B t ) denotes a Brownian Motion independent of ( W t ). Summarizing 


)• 


p (Cfc < iff(a)At, X tk > iff(a)At) < E ip(X t fc )l{x tjt>s(a)v ^} ' 
and we have (5.8) for Ak = {X tk > x(a)y/a). 

Step 2. Let us prove that for all x > x(a)pa: 

(1 - K (a)At) 2 (5(a)V^) 2(1 -“)' 


ip(x) < exp — 


2a 2 At 


(5.9) 


Applying again (5.5), we have 


= 1 erfc / [(1 — K( a )At)x + y/ab a (a)At) \ 

W 2 V V2Atax* ) 


exp 


x erfc 


2[bff(a) — K(a)x][x — (1 — y/a)b a (a)At] 


1 


Vv^A tax'* 
A(x) + B(x). 


[(1 + K{a)At)x — (2 — \fa)b a (ot)Ai\ 


26 
















Since At < l/(2A'(a)), and erfc(z) < exp(— z 2 ) for all z > 0 we have 


A(x ) < - exp I - 


[(1 — AT (a)At)x + y/ab <J (a)At\ \ 1 


< ^ exp ( - 


2o 2 Atx 2a 

On the other hand, for x > x{a)y/a, and At < l/(2A'(a)), it follows 

x > (2 — \/a)b cr {a)At /(1 + KAt) 
so the argument of the function erfc in B is positive, and then 


(1 — K{a)At) 2 x 2 ^ 1 
2a 2 At 


B(x) < i exp 


x exp j 
1 


2 \b a (a) — K(a)x\[x — (1 — \fa)b (T (a)At\ \ 


[(1 + K(a)At)x — (2 — y/a)b r7 (a)At\ 


= 2 eXP \- 


2Ato 2 x 2a 

[(1 — K(a)At)x + \Jab a (a)Atf 
2Ata 2 x 2a 


< 2 exp l 


(1 — K(a)At) 2 x 2 ( 1 


2<j 2 At 


ip(x) = A{x) + B{x) < exp ( — 


So 


and since the right-hand side is decreasing on x, we have (5.9). 

Step 3. Let us conclude. Putting together (5.8) and (5.9) we have 

N-l 

P(0 Q < T) < £ E (v-(X t Jl {3Ptfc>s(Q )^ } 


(1 — K (a) At) 2 x 2 ^ 1 a ) 


2a 2 At 


k =o 
N-l 

< ^ exp ( - 
fc=o 

C 


; }) 


(1 - K{a)At) 2 {x{a)^) 2 ^-^ 


2a 2 At 


’ (X tk > x(a)y/a) 


with 7 = (x(a)^/ r a) 2< ' 1 ~° , ' > /( 8 er 2 ). 


□ 


Proof of Lemma 2.4. From (1.4), standard arguments show that E[AS^(X) 4 ] < C{T). On the other hand, thanks to 
Corollary VI.1.9 on Revuz and Yor [18, p. 212], we have almost surely 


r T A0 a 


-kTAO Q PO — “ / 

e 4-U £ J o 


l M (I s )rf(I) 9 =0, 

because fore < (1 — i/a) 6 cr (a)Af, and s <T A 0 Q , l[o iE )(X fl ) = 0. a.s. Now, since 
L° T (X) = L^(X)t^Q a< T} + ^TA 0 Q (^)1{t<0„} = 


we can conclude that 


E [L° t (X) 2 ] = E[4(A) 2 1 {6q<t} ] < v /E[A°(A) 4 ]P( 0 a <T) < C^exp(-^). 


fl 
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5.4 On the negative moments of the stopped increment process (Z tA e j 

To prove Lemma 2.5 (see section 5.5 below), we need to control the negative moments of the stopped increment 
process {Z t ^e a }o<t<T- This is the object of the following lemmas, that can be summarize in the following 

Lemma 5.2. Let q > 1. Let © Q be the stopping time defined in (5.6). Let us assume At < A max (a). Moreover, let 
us assume 6(0) > 2ct(l — a) 2 when a £ (|, 1), and 6(0) > |er 2 ((7 + 1) when a = Then there exists a constant C 
depending on 6(0), a, a, T and q but not on At, such that 


vt e [o,T], 


E 



< C 



Existence of Negative Moments. Case a = \ 

The proof of the existence of Negative Moments of Z t A e a has two parts. First we study the quotient X v ^/Z s , and 
then we proof the main result of the section. 

Lemma 5.3. For a = \, and At < 1/(4 K) A Xq we have 


sup P 

0 <s<T 



< C At 8^2 


M1/2) 


(5.10) 


To prove this lemma, we need the following auxiliary result, the proof of which is postponed in Appendix A as a 
straightforward adaptation of the Lemma 3.6 in [5]. 

Lemma 5.4. Assume Hypothesis 1.1 holds, and 6(0) > cr 2 /4. Assume also that At < 1/(4A^) A Xq. Then, for any 
7 > 1 there exists a constant C depending on the parameters 6(0), K, a, Xq, T, and also on 7 , such that 


sup Eexp ( — 

k=0,...,N 


X 


tk 


OA , <C — 

yo A At) \ Xq 




Proof of Lemma 5.3. We start by proving 

X 


sup P I Z s < —] < sup E exp ( — X tk 


0 <s<T \ 2 / fc=o,. 1 ),Af ^ V 7 tj2 At J 

Indeed, if we call As = s — rj(s), and A W s = (W s — then 

plz s < ^ j < P ( a^X^AWs + 6 <J (l/2)As + (1 - KAa)X v{a) < 


< E 


AW S < 6 j (1/2)As + (4 -KAs)X v{b 
\J As 


tt \l X 

1 TV A \2' 


X, 


ri(s) 


< E exp — 


< Eexp — 


(frg(l/2)As + (| — KAs)X v ( s -)) 
2<j 2 AsX. 

(1 - 2KAt) 2 X v(s) 


v(s) 


8cr 2 A< J 

From here, the bound (5.11) follows easily, and then we conclude using Lemma 5.4. 


(5.11) 


□ 
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(5.12) 


Lemma 5.5. Let 0 1 be the stopping time defined in (5.6), and q > 1. If At < A max (l/2), and 

6(0) > ^ 2 (9+l)- 

Then there exists a constant C depending on 6(0), er, a, T and q but not on At, such that 

1 


Vi€[0,T], E 


Z. 


f A0 i 


<C 1 + 


Proof. Let us call A W s := (W s — W rl ( 9 f), and As := (s — t]{s)). By Ito’s formula 

r tAe h b(x v(e) ) 


But, 


Indeed, 


E 


7 q 
^tA 0i 


= \ ~ d E 


/ o Z? +1 


-ds 


+ ^ + 1) E 


ftA0 


$ 1 


/o z! +2 


( ' a \l X r,{s) + ds 


/ - - 2 \ 2 

( a \j X v ( s ) + — AWgJ < (t 2 A I) ( s ) + crZ s , P — a.s. 


(5.13) 


(5.14) 


— er 


\/+ “All's I — cr ( X v ( s ) a\/X rt ( s )AW s d—— AW t 


V(s) 


‘ill)' 


—o' + -A — (X vW + b(X v ^)As —— As) 

But, thanks to the Lipschitz property of b , 

2 2 
X-rj(s) + HX v (s))As —— As — ^(s) + (6(0) — KX V(S) ) As —— As 

= 6(t(1/2)As + (1 - KAa)X v(t) > 0, 

since As < At < 1/(2 K), and 6 CT ( 1/2) > 0. So we have (5.14). Introducing (5.14) in (5.13), and using b(x) > 
6(0) — Kx, we have 


E 


Z 


t A0 i 


— (7 - d E 


r tA0 


i 6(0) 


z! +1 


ds 


■ qK E 


[ tA@ i 

/o z ! +1 


ds 


q(q + 1) 


rtA@ i 

'“ E / —g+2 {^(«) + ■Z's} ds 

do Z„ 


(5.15) 


Since 


^ X v (s) ^ , r, 

Z, - Z s V.<*„e,/2}+ 2 ’ 


and applying Holder’s Inequality for some s > 0, we have 


E 


Z 


t A0 i 


<4t - # 


ftA0 


2 6 ( 0 ) 

zl +1 


ds 



' ,tA0i x 

+ 2gATE 

/ fff* 

do 


3q(q + 1) 


2 e 


ptA© 


2 1 

z! +1 


ds 


+ 


C 


A t?+ 2 


(EllJJ^P^a,^)/^ 1 £ ds. 
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Since b( 0) > 3 cr 2 (q + l)/2, we have 156 cr (l/2)/8tJ 2 > 2q + 2, so choosing e = q/(2q + 2), and applying Lemma 
5.3 we have 

P(Z s <X, (s) /2) 1 - £ <CA^+ 2 , 

and then 


E 


7~1 
J tA@ 1 


1 


2qKE 


<-t A 0 i 




G?S 


/ 3(g + 1) 2 


ff 2 - 6(0) E 


rtAQ i 


z ! +1 


ds 


C. 


Since from the Hypotheses, the third term in the right-hand side is negative, we can conclude thanks to Gronwall’s 
Lemma. □ 


Existence of Negative Moments. Case a > \ 

Lemma 5.6. For «£ (), 1), ifb( 0) > 2a(l — a) 2 er 2 and At < l/(4aA'(a)), there exists 7 > 0 such that 

„<”?/ ( z - 5 ( x “ g) x -«) s exp (“St) ■ f5 - 16) 

Proof. Let us call A W s := {W s — 1^^)), As := ( s — rj(s)), and 

q(X v{s) ,AW s ) = ^X^AW 2 + aX; {s) AW s + ^ + (M«) - K(a)X v(s) ) As. 

Notice that for fix x £ K, q(x, •) is a quadratic polynomial. Using (5.2), we have 

p (Zs < ^1 - ^(s)) < P(q(Xr,( s ),AW s ) < 0,x v(s) < 5(a)) 

+ p{q(X vis) ,AW s ) < 0 ,x v(s) > 5 (a)), 

where recall, 5(a) = b a (a)/K(a). But 


q( x v (s),&W s ) < 0,X V{S) < x(a) 


= E 


P (q(x, VAsJ\f) < 0 


=^ (a) ' IL { x 'iW^ s(a )} 


where Af stands for a standard Gaussian random variable. As in the Lemma 2.3, we have a quadratic polynomial in 
Af, its discriminant is 

A = a 2 x 2a As — 2aer 2 a ; 2a ^ 1 As + ( b a (a ) — K(a)x) As = —2aa 2 x 2a - 1 As 2 (b c r (a) — K(a)x ), 

. 2 a 

so if x < 5(a), A < 0 and the quadratic form in Af has not real roots, and in particular is non negative almost surely. 
Then 

P(q(X v{s) ,AW s ) < 0,X v{s) < 5(a)) = 0. 

On the other hand, 

v(q{ x v (s),AW s ) < 0, AC t?(s) > 5(a)) 


< E 


PLV<- 


b a (a)As + (X - K(a)As) x 


c=X 


■n(s) 


ax a V As 

and since At < 1/(4aA'(a)) we can apply the exponential bound for Gaussian tails and get 

“ K{a)As) 2 x 2 ^ 1 -^ 


{^r,(s)>a:(a)} 


p(?(^r,( s ), AWs) < 0,A, (S) > a;(a)) < E 
We conclude by taking 7 = 5(a) 2 ( 1_Q )/(16er 2 ). 


exp -- 


er 2 As 


c=X 


\-W( s ) — x ( a )} 


rj(s) 
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Lemma 5.7. Let 0 a be the stopping time defined in (5.6). Let us assume for a 6 1), 6(0) > 2a(l — a) 2 , and 

At < A max (a), then for all q > 1, there exists a constant C depending on 6(0), a, a, T and p but not on At, such 
that 


Vt £ [0,T], 


E 



< C 



Proof Let us call A W s := W s — W v ( s y By Ito’s formula and the Lipschitz property of 6, 


E 


Z. 


t a0 q 


< ^ - qE 


ftA0 Q 


m 

zl +1 


ds 


+ qK E 


ft A0 q 


A 


rj(s) 


zl +1 


ds 


g(g + 1) 


E 


/•t A© q 


1 / 


x! +2 




(5.17) 


Following the same ideas to prove (5.14), for all s £ [0, i] we can easily prove that almost surely 

(aX; {s) + aa 2 xl7^AW s Y < +X 2 “, + 2a+X^X,. 


L v( s ) 


vi s ) 


L r;(s) “s- 


Introducing this bound in the previous inequality, we have 
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Z 


tA0 o 


<~q ~ # 


f tA0 Q 


6 ( 0 ) 

X! +1 


ds 


■ gXE 


ftA0 Q 


X 


7)(s) 


x! +1 


ds 


+ *±il^ E 


rtA@ a 


1 f^2a 


/o X? +2 


f cr -rr—2a— 1-^ 1 . 

|-^7j(s) + 2 q ++(s) Z S J « s 


(5.18) 


since for r £ {1, 2a — 1, 2a}, 


X 


T)(S) 


< 


X 


rj(s) 


2 a 

~{z a <x n(s) ( i-|a)} + l, 2a - 1 


we get 


E 



1 

[ f ^ 6(0) dsl 

[Jo zr \ 

2 a 

ftA0 a 

1 

+ tA0 a J 

<-, - gE 
^0 

+ 2a —1^ E 

1/ 



g(g + 1 ) j2 (2a) 2a+1 E 


2 ’ (2a - l) 2a 

tA0 a 


ft A0 q 


-g+2(l-a) 


ds 


CE 


/o 


Xri(s) 

x! +1 


—2a -r^a-l ' 

x! +2 x! +1 


1- _ i „ds 


The last term in the previous inequality is bounded because of the definition of 0 Q and the Lemma 5.6. Indeed, 
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ftAOo 
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7j(s) 
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2 a 
V («) 


x; 


2a—1 ■ 

v( s ) 


?q +1 
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c ‘- T 

A f?+ 2 
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9+1 


1 {Z a <^ ( , ) (l-ia)} dS 


Af<?+ 2 


r 

exp (ii) 


v (») 


“r 2a i t 720-1 
A lW +A r(W 


) 2 ]p(z,<x, w (i-T)) 
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ds 


31 



























































So, (5.18) becomes 


E 


Z. 


tA0 Q 


<—q - # 


rt a© q 


/o 


6 ( 0 ) 

z ! +1 


ds 


2a 


d(<7 + 1) _ 2 (2a) 2 “ +1 
+ CT ( 2 a-l) 2 “ E 

Ai A 2 


2a - 1 
ftA0„ 


5 ATE 
1 


rt AO a ^ 

/ 

/0 A, 


/o z 9+2(1 ” Q) 


ds 


+ C. 




< —w + 2qI\E 

rtAOc, ^ 

E 

7 9 

^tAOc 

/ 



x o 

Jo z s 


But, for any Ai, A 2 > 0, the mapping 2 1 —> _ q+ ff_ a) — is bounded, and (5.19) becomes 

ptAOc 

■C, 

y° z s 

from where we can conclude applying Gronwall’s Lemma. 

5.5 On the corrected local error process 

Proof of Lemma 2.5. Let us recall the notation in the proof of the main Theorem 

D S (X) := ai: - aX a v{s) - ao*X^\w a - W v{s) \ 


and also introduce 


— —ct 2 —2 ct —1 

SuaoA*) -vXni sAO a ) + aa X V (sA& a )O^uA0 a - W„(.Ae o) ), 


and AW S := (W a - W„ (a) ). 

Using Lemma 5.1, and the hniteness of the moments of I), is easy to prove 

E [D s (X) 2p ] < CE [D sA e a (X) 2p l {ea > v{s)} ] + CAt 2p . 

Then we only have to prove 

E[D sA e a (X) 2p he a >r,( s )}]<CAt 2p . 

Notice that A' sA o q = Z sA Q a , so 

— I —ol — ot. 2 — 2(x — 1 | 

t9 S A0 c ,(A')l{© Q ,> 7) ( s )} = | crZ sA Q a - o-A^ ( sA0q ) - acr A’„( aAea )AWg A e a | l{0 a >??(s)}- 

Then applying Ito’s Formula to the function <r\x\ a which is C 2 for x > CAt. we have 


-DsA0 q (-A)1{©„>77(s)} = 


sA0 o 


a ^ri(s AO a )dW u 


1 ~ry^~ a a 

v(sAS a ) \Z uA Q a A r/(sA0Q 
rsA0 a rv 2 rr 3 ~F 2a — 1 

+ / - _J 1 ^ S A 0 Q ) A w uA0a dW u 

Jv(sAO a ) Z uAQot 

PSAQ a 

I _ 1 —a K-^i?(sA@ q ))^ 

Jn( ' ■ " 


'v(sAO a ) Z uA Q a 


pSAQr, 


1 a(l — a)cr 


/ 9 ^2 — 0 ; 
,r l(sAQ a ) Z Z uAQq 


SuAe a (X) du j* 1 {© q ,>?7(s)} 


—• J\ J 2 H” *^3 — ^4- 

Notice that on the event { 77 ( 5 ) < 0 Q } we have 77 ( 5 ) = 7y(s A ©a), and then 

nSAQ a / _ _ ^ ^ 

K{») 


E[|J!| 2p ] =E 


1 (e Q >r,( S )} _i_ a - =13 


^UAOa Xrj(sAOa) , 


Cr '^r;(sA 0 a ) 


2p' 


(5.19) 
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(5.20) 


(5.21) 


(5.22) 
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By the Burkholder-Davis-Gundy inequality, there exists a constant C p depending only on p such that 
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aa 


7 * “ y 1 a 
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—i-Q ^rl-a 
^«A0 a ^ rj(sA<d a ) 


observing that the integrand in the right-hand side is positive. And we have 

E[|Ji| 2p ] < (aa 2 ) 2p C p E 


4a _ 2 


< CE 


v(s) 
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?|(sA0 a ) J1 {0a>>)(s)} 
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^(sA0 a )^{0a >77(5)} 


But for x, y > 0, and /? € [0, it holds |a:^ — y l3 \(x 1 & + y 1 < 2\x — y\, so 

E[|Ji| 2p ] < CE 


/ r® X^ a 2 "a p 

I I rv a 2 n A »)(s A e a ) 

I / ^^-r 7 (sA@Q.) ^uA© a ) *-{©a>v( s )} _ 2 

Z„ A e o 


du 


E 




< CAt?- 1 

Let a > 1. Thanks to Holder’s inequality we have 

,2p 


-=2p(2a-l) 

(-y -7 \ 2 P t A p(s A 0 a ) 
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7 2 P 
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~y 2a P 

^«A 0 O 
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We use Lemma 2.1 to bound the Local Error of the scheme 
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(^ij(s A e 0 ) — Z uA Q a ) ( “ 11 Iq©^ 


»)(«)} 


< CAf^r, 


On the other hand, when a > 5 , we have control of any negative moment of so 
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<C, 


whereas when a = |, we choose a > 1, such that 26(0)/cr 2 > 3(2 ap + 1), so we have control of the 2ap-th negative 
moment of Z uA Q a . And then 

—Aap(2cx—1) 

' = E 
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So, in any case we have 
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(^77(sA0c) ^U/\Qa) ^{O a >?7(s)} 


—2p(2a-l) 
i)(sA0„) 

-y 2 P 


< CAt p . 


And then we can conclude E[| Ji| 2p ] < CAt 2p . 

Using the same arguments for E[| J 2 \ 2p ], we have 


K[\J 2 \ 2p } < C P E 


< CAt p ~ 


< CAt p ~ 
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/t)(s) 


lri(s] 


w? 2 - 1 ’ 

zi 

2(2a-l)p 


rsAO a a 2 a 6 XZ) a 2 
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l{e a >, W} ^A "fte* 
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>7(s) 

—4(l-a)p 
^uAQc 


To bound E[| J 3 | 2p ] we proceed as follows 


E[|J 3 | 2p ] = e 


f sAea aa — ' 

/ ^{0 O >'/W}—iTi n(sAO a ))d’U 

’«') Z U Ae a ) 


® ( ^{ 0 c>??(s)} AW% @1 )du<CAt 2p . 


2 P 


< (acr) 2 p Af 2p_1 f E (—^ 6 (^ ( . A e o) ) 

\ 2 «A0 q 
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du 


1 —a 


< {aa) 2p At 2p ~ 


< CAt 2p . 
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/?j(s) 


XuA0 a 


E [ b ( X r,(s A0«l) ” 


Finally for E[| J 4 | 2p ] we consider first a > A . In this case we have control of any negative moment of Z uA ® a . So 


proceeding as before 


E[|J 4 | 2p ] =E 


(*sA0 q 


lrt(s) 


l{0 a >r)(s)} 2 — 2 - 
Z uAQ 


la(l-a)a V2 , \ 

7 t — 2 -Q <Sii A 0 a (X) du ] 


\ 2 p 


< CAi 2 '- 1 Z E a* 


< CAf 2p . 

The case a = \ is a little more delicate. Let us recall the identity used in the proof of (5.14) 


vV2 




^,(M0i) + yAW^ABiJ =o-Z„ Aei -(T [b{X v{sAe Z~ — J (uAQi -7)(sA0i)), 

so, we have from the definition of 01 
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SuAsZX^ = (aXZ A&i ) + ^AWuaq^ < c (zt P A@i + (b(X v(sAe ^ ^ ? At 2 ^ 

<C^l+^(X, (sAei) )-y) ^ 


2 V 

uAQ i 
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Then 


E[| J 4 | 2p ] < CAt 2p ~ 1 

1 E 


1 v{ s ) 

< CAf 2p_1 

[ E 


d r]( s ) 

< CAf 2p_1 

f \ 


J V (s) \ 

< CAt 2p . 



-=3p ( ff ^i)(sA0i) 
. Z uAe a v 2 


—AW uA ei 

Z 2 


4p 


dw 


; uA0 q 


1 + ^b(X v ( s A Qi)) 



E 




^(^i)(sA0i)) 



So for every a€ [^, 1), E[</ 2:p ] < CAt 2p , from where we conclude on the Lemma. 
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A Appendix 


Proof of Lemma 1.4. Let us recall the notations As = s — rj(s), and A W s = W s — I-F,^.,). Let us define r m = inf {t > 
0 : X t > m}. Then by Ito’s Formula, Young’s inequality and the Lipschitz property of b, we have 

‘ ft/\Trr, 


E[X t 2 ^J < xg + CE 


2 P 


X* p + C + X% s) ds 


+ CE 


rt/\T m 2^2ct—1 . \ 

{ aX V (s) + acr X r,(s) XW s) ds 


From the definition of X, a straightforward computation shows that for all s £ [0, t] almost surely 


x] p < C ( 1 + X% s) + AWf~ a + (A W 2 - As) 

Putting this in (A.l), we have 


E[^(Ar ra ] < 4 P + CE 


CE 


J 1 + X n(s) + {fi jX r t (s ) + a<jZX ri(s) AW S ) ds 
™ A + ( AW 2 - (s - r?(s))) ds 

' ptA Tm 

f 1 + X T( S ) + X ?(s) + x nt) AW 2p ds 

.JO 

(AW 2 - As) ds . 


2 v 


< x 2p + CE 


Since a £ [^, 1) we have < 1 + X^ s ), and then, using Young’s Inequality and the finiteness of the moments 

of Gaussian random variables, we conclude 
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/ E 
Jo 

AW, 1 - 

+ E 


^2 p 


(A.l) 


n X tArJ < CX 2P + CE 
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<Cx 2p + C su P E[X* AT Jds. 
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Since the right-hand side is increasing, we can take supremum in the left-hand side and from here, applying Gronwall’s 
inequality, and taking m —> oo we get 

supE[Jf( P ] < Cx g P . 

t<T 

From here, following standard argument using Burkholder-Davis-Gundy inequality we can conclude on Lemma 1.4. 

□ 


Proof of lemma 5.4. First, from the definition of X tk we have 

X tk > X tk _, + (M 1/2) - KX tk _ 1 )At + (W tk - W tk _ k ), 


then 


Eexp(-/u 0 X t J <Eexp(-/^ 0 X^, + {b a (\/2) - KX tk _f)At + X tk _ x (W tk - W tk _ t ) 
where = l/'yo 2 At. From here, just as in Lemma 3.6 in [5], we conclude 


Eexp (~fj. 0 X tk ) < exp (—/r 0 6 cr (l/2)Af) Eexp ( ~n 0 X tk _ 1 




1 - KAt - 


a 2 At 


-Po 


(A.2) 


Then if we introduce the same sequence (pj)j > 0 of Lemma 3.6 in [5], given by 


Pi = 


'3)3 
1 

7<t 2 A t 5 

Pi -1 


1 - KAt - ^Pj -1 


3 = 0, 

3> 1- 


We can repeat the proof in [5] and find out that if At < 1/(2K) then, the sequence (pj)j > 0 is nonnegative, 
decreasing and satisfies the following bound 


Pi > Pi 


1 


1 + ^-A t(j - l)/i 0 


- K 


Xt{j - l)/r 0 
1 + ^-At(j - l)/r 0 


, Vj > 1. 


On the other hand making the same calculations to obtain (A.2) we can get for any j £ {0, .... k 1}, 


Eexp (- PjX tk _.) < exp (-n j b rT (l/2)(\)At) Eexp ( -p J X tk _._ 1 

from where, by an induction argument we have 


1 - KAt - 


a 2 At 


-Pi+i 


k -1 


E (—poX tk ) < exp ( -b a ( 1/2) E PjAt ] exp (x 0 /Ltfc). 

i =o 


From here, and the bound for the sequence (// ; ) ;/ > 0, we have 


E(-«,X„)<C( 


2b a (l/2) 




From where we see immediately 


sup E exp- 2 ~aT 

k=o,...,N V 70 " A t 


<C[ — 
x o J 


A t\ ^ ^ 27; 


. □ 
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